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Abstract 

A  further  analysis  of  the  IFRF  semi-industrial- scale  experiments  on  flameless  (mild)  combustion  of  natural  gas 
is  carried  out.  The  experimental  burner  features  a  strong  oxidizer  jet  and  two  weak  natural  gas  jets.  Numerous 
publications  have  shown  the  inability  of  various  RANS -based  mathematical  models  to  predict  the  structure  of  the 
weak  jet.  We  have  proven  that  the  failure  is  in  error  predictions  of  the  entrainment  and  therefore  is  not  related  to 
any  chemistry  submodels,  as  has  been  postulated. 

©  2007  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

Our  paper  is  concerned  with  rapidly  develop¬ 
ing  combustion  technology  [1]  that  has  been  given 
different  names;  high-temperature  air  combustion 
(HTAC  [2]),  mild  combustion  [3-5],  or  flameless  ox¬ 
idation  [6].  The  essence  of  this  technology  is  that  fuel 
is  oxidized  in  an  environment  that  contains  a  substan¬ 
tial  amount  of  inert  (flue)  gases  and  some,  typically 
not  more  than  3-5%,  oxygen.  For  clean  gaseous  fu¬ 
els  that  do  not  contain  any  fuel-bound  nitrogen,  this 
results  in  very  low  NOx  emissions  even  if  the  com¬ 
bustion  air  stream  is  preheated  to  temperatures  in 
excess  of  1000  °C. 

There  are  several  methods  of  realization  of  this 
combustion  technology,  as  underlined  in  Refs.  [1,4], 
and  two,  perhaps  the  most  common  ones,  are  de¬ 
scribed  below.  In  the  first  method  (see  Refs.  [1,2]), 
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the  combustion  air  is  provided  by  a  central,  strong 
(high-momentum)  air  jet  that  is  surrounded  by  a  num¬ 
ber  of  weak  (low-momentum)  fuel  jets  (in  industrial 
applications  typically  two  jets  are  used).  These  weak 
fuel  jets  are  positioned  away  from  the  central  air  jet  in 
order  to  inject  the  fuel  into  recirculated  combustion 
products  and,  by  doing  so,  to  dilute  the  fuel  before 
it  mixes  with  the  combustion  air  stream.  The  second 
method  is  recognized  in  the  literature  as  a  “classical” 
method  of  achieving  flameless  combustion  [6].  It  can 
be  characterized  by  a  central  fuel  jet  and  a  number  of 
air  jets  positioned  in  the  relative  vicinity  of  the  central 
fuel  jet. 

Predicting  performance  of  industrial  devices  that 
utilize  the  mild  combustion  principle  is  complex, 
since  in  addition  to  considerations  on  combustion 
chemistry,  the  fluid  flow  and  mixing  between  fuel, 
oxidizer,  and  combustion  products  have  to  be  consid¬ 
ered.  The  publications  leading  to  both  the  develop¬ 
ment  and  validation  of  such  models  can  be  split  into 
these  concerned  with  the  strong  jet-weak  jet  design 
and  those  concerned  with  a  multiplicity  of  jets  of  sim¬ 
ilar  momentum.  In  conjunction  with  the  work  of  the 
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first  group,  a  unique  set  of  measured  data  was  gen¬ 
erated  [1,7].  The  availability  of  these  data  initiated 
a  number  of  publications  aiming  at  both  the  devel¬ 
opment  and  validation  of  several  mathematical  mod¬ 
els.  In  all  these  publications,  the  Reynolds-averaged 
Navier-Stokes  (RANS)  approach  was  retained;  how¬ 
ever,  the  works  differed  in  the  turbulence  models  or  in 
turbulent-combustion  interaction  submodels  used  or 
in  both.  In  Mancini  et  al.  [8],  the  standard  k-s  turbu¬ 
lence  model  was  used  together  with  three  combustion 
submodels:  the  eddy-breakup  model  with  a  two-step 
reaction  scheme,  the  eddy-dissipation  concept  model 
with  chemical  equilibrium,  and  the  pdf/mixture  frac¬ 
tion  model  with  nonadiabatic  lookup  tables  (for  de¬ 
tails  of  the  NO  postprocessor  see  Ref.  [8]).  With  the 
exception  of  the  fuel  jet,  the  computations  resulted  in 
predictions  of  good  quality.  The  predicted  NO*  fur¬ 
nace  emissions  were  in  very  good  agreement  with 
the  measured  values  for  a  number  of  combustion- 
model/NO*  postprocessor  combinations.  However, 
the  models  failed  to  predict  the  structure  of  the  weak 
jet.  More  precisely,  within  the  fuel  jet,  the  computed 
temperatures  were  typically  500-600  °C  lower  than 
the  measured  values  and  the  predicted  CH4,  CO2,  and 
CO  concentrations  differed  by  a  factor  of  2-3  from 
the  measured  values.  Dong  [9]  used  a  RANS  model 
based  on  k-s  equipped  with  a  mixture  fraction/pdf 
submodel  to  simulate  the  same  experiments  [1,7]  and 
his  predictions  suffered  exactly  the  same  deficiency. 
More  recently,  Awosope  et  al.  [10,11]  used  a  Borghi 
diagram  [12]  to  identify  the  combustion  regimes  of 
the  considered  experiments  [1,7].  They  used  a  com¬ 
bustion  submodel  that  included  the  transport  equa¬ 
tions  for  the  mean  mixture  fraction,  its  variance,  and 
the  mean  reaction  progress  variable.  The  source  term 
in  the  reaction  progress  variable  was  modeled  follow¬ 
ing  the  work  of  Lindstedt  and  Vaos  [13],  which  in¬ 
cludes  the  influence  of  strain  on  the  flamelets  by  small 
eddies.  The  flamelet  functions  were  obtained  from  a 
priori  computed  lookup  tables.  Similarly  to  the  mod¬ 
els  of  Mancini  et  al.  [8]  and  Dong  [9],  the  model  of 
Awosope  et  al.  [10,1 1]  also  failed  to  predict  the  struc¬ 
ture  of  the  weak  fuel  jet.  Notwithstanding  the  fact  that 
perhaps  a  more  thorough  analysis  of  the  applicability 
of  the  above-listed  chemistry-turbulence  models  of 
the  IFRF  experiments  [1,7]  is  required,  as  initiated  by 
Awosope  [10,11],  one  cannot  refrain  from  observing 
that  the  more  sophisticated  the  chemistry-turbulence 
model  is,  the  worse  are  the  predictions  of  the  low- 
momentum  fuel  jet. 

2.  Objectives 

The  goal  of  this  paper  is  to  resolve  the  issue  of 
the  poor-quality  predictions  for  the  weak  fuel  jet. 


The  main  question  is  whether  (new)  low-temperature 
chemical  kinetics  schemas  are  really  needed  to  im¬ 
prove  the  predictions  in  this  flame  region.  Such  a  sig¬ 
nificant  discrepancy  between  the  predictions  and  the 
measurements  in  the  fuel  jet  should  have  resulted  in 
error  predictions  of  the  furnace  exit  NO  emissions. 
To  the  contrary,  all  the  NO  postprocessors  have  pre¬ 
dicted  the  furnace  NO  exit  values  with  relatively  good 
accuracy.  To  explain  this  rather  surprising  mathemat¬ 
ical  model  effect,  we  use  the  IFRF  measurements 
[1,7],  several  RANS  models,  and  the  radical  chem¬ 
istry  mechanism  (GRI  Mech  3.0  [14])  incorporated 
into  the  Chemkin  code. 


3.  The  measured  data 

The  in-furnace  measured  data  used  in  this  paper 
originate  from  the  IFRF  experiments  [1,7]  carried  out 
in  a  refractory  lined  furnace  with  a  2  x  2  m  cross  sec¬ 
tion  and  a  length  of  6.25  m  (see  Fig.  1  in  Ref.  [7]). 
The  furnace  was  equipped  with  one  (0.58-MW  ther¬ 
mal  input)  burner  operated  under  steady  state  condi¬ 
tions.  The  burner  consisted  of  a  central  oxidizer  jet 
(named  also  comburent  jet)  and  two  natural  gas  in¬ 
jectors  positioned  28  cm  away  from  the  central  jet. 
A  gaseous  stream,  of  temperature  1573  K,  contain¬ 
ing  19.5%  wet  O2,  59.1%  wet  N2, 15%  wet  H2O,  and 
6.4%  wet  CO2,  was  supplied  into  the  furnace  through 
the  central  jet.  This  oxidizer  stream  (830  kg/h)  con¬ 
tained  also  94  vppm  of  NO  calculated  on  a  wet  basis 
(110  vppm,  dry).  The  volumetric  composition  of  the 
furnace  exit  gases  was  1.6%  O2,  27.3%  H2O,  13.7% 
C02,  57.7%  N2  and  102  ppm  of  NO  (140  vppm, 
dry).  The  “net”  NO  production  rate  in  the  furnace  is 
0.06  g/kg  of  fuel.  Detailed  mapping  of  the  turbulent 
velocities,  temperatures,  and  chemical  composition 
(02,  CO,  C02,  CH4,  H2,  NO,  and  NO*)  as  well  as 
radiation  has  been  reported  elsewhere  [1,7]. 

4.  Mathematical  models 

4.1.  The  RANS  mathematical  models 

For  predicting  the  measured  data  we  use  several 
RANS  models.  The  turbulence  is  modeled  using  ei¬ 
ther  a  standard  k-s  model  or  a  RNG  k-s  model 
(more  elaborate  Reynolds  transport  models  were  also 
tested  and  only  minor  changes  were  observed).  We 
use  three  combustion  submodels  listed  in  Table  1. 
The  eddy-breakup  submodel  with  a  two-step  reaction 
scheme  [15]  is  the  most  simple  one  (model  1  in  Ta¬ 
ble  1).  The  second  submodel  (the  standard  Fluent  5 
CFD  code)  is  based  on  a  mixture-fraction  transport 
equation  that  is  equipped  with  nonadiabatic  lookup 
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Table  1 


Measured  and  predicted  values 


Measurements  [1,16] 

Predictions  using  model  no./run  no. 

12  3 

4 

5 

Furnace  exit  temperature  (K) 

1493 

1519 

1498 

1489 

1480 

1489 

Furnace  exit  O2  (%,  vdry) 

2.2 

1.8 

2.14 

2.2 

2.1 

2.19 

Furnace  exit  NO  (ppmvd) 

140 

137 

130 

129 

140 

135 

Furnace  exit  CO  (ppmvd) 

0 

0 

8 

10 

15 

13 

Peak  temperature  (K) 

1808 

1948 

1973 

1922 

1770 

1840 

Peak  H2  (%,  vdry) 

1.97 

9 

13.1 

10.9 

2.5 

2.0 

Peak  CO  (%,  vdry) 

3.23 

21 

33 

28 

3.1 

2.8 

Note.  1 — IFRF-EBU  [15]  model;  2 — Fluent  5-pdf  (mixture  fraction)  [17];  3 — Fluent  5-pdf  (mixture  fraction) — fluctuations  of 
the  mixture  fraction  have  been  switched  off;  4 — reactor  network  model  with  GRI  3.0  mechanism  and  entrainment  calculated 
using  RANS  model- 1;  5 — reactor  network  model  with  GRI  3.0  mechanism  and  entrainment  estimated  from  the  measured  data. 


Natural  Gas 
Jet 


R  22 


Fig.  1.  Series  of  perfectly  stirred  reactors  representing  the  strong  jet  (reactors  R1-R10)  and  the  weak  jet  (reactors  R12-R21). 


tables  (model  2  in  Table  1).  In  the  third  submodel, 
the  turbulent  fluctuations  are  switched  off,  but  other¬ 
wise  its  formulation  is  identical  to  model  2.  To  each 
of  the  above  combustion  submodels  the  same  NOx 
postprocessor  is  appended.  The  thermal  and  prompt 
NO  are  calculated  in  the  same  way  as  in  Ref.  [15], 
with  no  changes  to  any  of  the  model  parameters  and 
constants.  However,  an  extra  NO-reburning  path  fol¬ 
lowing  the  work  of  Chen  [16]  and  De  Soete  [17]  was 
added,  as  explained  in  Ref.  [8].  Although  we  use  the 
Chen  NO  reburning  model,  we  realize  that  in  the  light 
of  the  present  knowledge  of  NO  precursors  [18]  the 
model  does  not  have  a  firm  scientific  basis. 

4.2.  Series  of  perfectly  stirred  reactors 

To  provide  the  background  for  the  detailed  analy¬ 
sis  of  the  chemistry,  we  have  subdivided  the  flow  into 
a  network  of  perfectly  stirred  reactors  (PSRN).  We 
have  divided  the  strong  central  jet  into  10  reactors 
marked  in  Fig.  1  as  R1-R10.  Each  reactor  is  of  the 
same  length  and  its  width  is  determined  by  the  po¬ 
sition  of  the  zero  axial  velocity,  which  is  determined 


by  the  jet  expansion  angle.  The  location  of  this  zero- 
axial-velocity  line  is  known  from  both  the  measure¬ 
ments  and  the  RANS  predictions.  The  weak  natural 
gas  jet  is  divided  into  10  PSRs  marked  in  Fig.  1  as 
R12-R21.  The  streams  leaving  the  fuel  jet  reactor 
R21  and  the  center  jet  reactor  RIO  merge  into  reactor 
Rll.  The  rest  of  the  furnace  is  represented  by  reac¬ 
tor  R22.  As  has  been  shown  in  Fig.  1,  reactors  entrain 
hot  recirculating  combustion  products  that  originate 
from  the  last  reactor  (R22).  The  inlet  to  the  R12  reac¬ 
tor  is  a  47/2  =  23.5  kg/h  flow  of  natural  gas  of  the 
following  volumetric  composition:  0.88  CH4,  0.05 
C2H6,  0.02  C3H8,  0.04  N2.  The  inlet  to  R1  reactor 
is  a  830/2  =  415  kg/h  comburent  flow  of  the  vol¬ 
umetric  composition  specified  previously.  We  stress 
here  again  that  this  stream  contains  94  vppm  (wet) 
NO.  The  network  of  the  reactors  is  solved  using  the 
CHEMKIN  3.7  computer  code  [19]  with  the  GRI  3.0 
reaction  scheme  for  methane/air  flames  [14]. 

As  shown  in  Fig.  1,  both  jets  entrain  combus¬ 
tion  products  of  reactor  R22  and  we  use  two  meth¬ 
ods  for  determination  of  this  entrainment.  Since  a 
RANS  model  solution  is  obtained,  a  simple  integra- 
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Fig.  2.  Entrainment  of  the  jets;  left — the  weak  natural  gas  jet,  right — the  strong  comburent  jet;  (•)  entrainment  estimated  using 
the  measured  data  [1,16] — the  error  bars  indicate  the  inaccuracies  in  estimating  the  entrainment;  (□)  entrainment  calculated 
using  RANS  model- 1. 


tion  provides  the  entrainment  values.  Fig.  2  shows  the 
entrainment,  for  both  the  weak  and  the  strong  jets, 
obtained  using  the  RNG  k-e  based  RANS  model- 1 
of  Table  1.  In  the  second  method  we  use  measured 
data  that  include  the  time-mean  and  fluctuating  veloc¬ 
ities  (laser  Doppler  anemometry)  and  the  time-mean 
temperatures  and  gas  composition.  These  data  have 
been  collected  point  by  point  at  several  traverses  in 
the  symmetry  plane  cutting  trough  the  middle  of  the 
furnace.  Using  these  data,  we  have  calculated  the  en¬ 
trainment  rates  for  both  jets.  Since  the  strong  cen¬ 
tral  jet  remains  axisymmetric  and  “stiff,”  with  well- 
defined  boundaries,  this  procedure  is  accurate  and,  as 
is  shown  in  Fig.  2,  the  entrainment  rates  calculated 
from  the  measured  data  agree  well  with  the  RANS 
predictions.  The  same  procedure  has  been  applied  to 
the  weak  jet  of  natural  gas.  Since  the  weak  jet  is 
strongly  tilted  toward  the  center  jet,  our  procedure  for 
calculating  the  entrainment  using  the  measured  data  is 
in  this  case  less  accurate.  The  accuracy  with  which  we 
estimate  the  entrainment  rates  is  shown  in  Fig.  2  using 
the  error  bars.  The  same  figure  shows  that  there  is  a 
substantial  discrepancy  between  the  RANS  model- 1 
predicted  (other  RANS  models  provide  similar  en¬ 
trainment  rates)  and  the  “measured”  entrainment. 

5.  Results  and  discussion 

Table  1  shows  the  comparison  between  the  mea¬ 
surements  and  predictions  in  terms  of  the  furnace  exit 
values.  Excellent  agreement  is  obtained  for  all  the 
models  used.  The  overall  mass  and  energy  balance 
(including  radiation  balance)  is  satisfied  with  an  accu¬ 
racy  better  than  0.5%  for  all  the  computational  runs.  It 
is  worth  noting  that  simple  chemical  equilibrium  cal¬ 
culations,  carried  out  for  the  furnace  inlet  conditions 
and  for  the  furnace  exit  temperature  of  1493  K  and 
1  bar  pressure,  result  in  NO  emissions  of  392  ppmvd, 
while  140  ppmvd  was  measured. 


Fig.  3  shows  the  predicted  temperature  and  oxygen 
in  both  jets  using  the  RANS  model- 1  and  the  reactor 
network  model,  while  Fig.  4  shows  the  NO  and  NO2 
concentrations  within  both  jets.  The  predictions  using 
RANS  models  2  and  3  (see  Table  1),  show  similar  re¬ 
sults  in  the  fuel  jet  and  are  not  plotted.  For  the  reactor 
network  model  the  entrainment  figures  are  essential 
model  inputs  and  therefore  two  computational  runs 
are  performed.  In  run  4  (Table  1)  the  entrainment  is 
obtained  using  the  RANS  model- 1  predictions,  while 
in  run  5  (Table  1)  the  entrainment  is  estimated  di¬ 
rectly  from  the  measured  data  (see  Fig.  2).  Thus,  in 
run  5  (Table  1),  the  entrainment  into  the  fuel  jet,  at 
the  point  of  merging  into  the  comburent  jet,  is  larger 
than  in  run  4  by  a  factor  of  3. 

The  central  jet  is  easy  to  predict,  as  shown  in 
Figs.  3  and  4.  Until  a  distance  of  around  70  cm,  which 
corresponds  to  the  Rll  reactor,  there  are  no  com¬ 
bustion  reactions  taking  place  in  the  jet  and  its  tem¬ 
perature  is  determined  by  the  energy  balance,  which 
includes  the  entrainment  and  the  radiation  losses.  The 
oxygen  content  in  the  jet  decreases  linearly  with  the 
distance  due  to  the  entrainment.  Similarly,  the  NO 
concentration  increases  slightly  with  the  distance,  and 
this  effect  is  again  exclusively  due  to  the  entrain¬ 
ment. 

Consider  now  the  weak  jet.  Fig.  3  shows  the  pre¬ 
dictions  in  the  fuel  jet  obtained  using  the  reactor  net¬ 
work  models  (runs  4  and  5,  Table  1).  The  flow  rate 
of  CH4  leaving  the  fuel  injector  is  5.055  g/s.  The 
CH4  flow  rate  at  the  exit  of  reactor  21  is  calculated 
to  be  5.0335  g/s  in  run  4  while,  for  run  5,  a  figure  of 
4.84  g/s  is  applicable.  The  oxygen  concentration  in 
reactor  21  is  calculated  to  be  1.21  and  1.3  vol%  for 
runs  4  and  5,  respectively.  In  run  4,  only  0.43%  of 
CH4  is  converted,  in  the  fuel  jet,  into  the  intermedi¬ 
ate  radicals  and  CO  (0.2%  CO),  while,  for  run  5,  the 
corresponding  figures  are  4.3%  of  CH4  and  0.933% 
CO,  respectively.  Thus,  there  is  hardly  any  “consump¬ 
tion  of  fuel”  taking  place  within  the  fuel  jet.  The 
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Fig.  3.  Temperature  (left)  and  O2  (right)  in  the  jets:  (•)  measurements  [1,16];  (O)  RANS  model-1  predictions;  (V)  reactor 
network  model  with  RANS  model- 1  computed  entrainment;  (□)  reactor  network  model  with  entrainment  estimated  from  the 
measured  data  (see  Fig.  2). 


Distance  (cm) 

Fig.  4.  NO  concentrations  in  the  jets:  (•)  NO*  measure¬ 
ments  [1,16];  (O)  RANS  model-1  predictions;  (V)  reactor 
network  model  with  RANS  mode-1  computed  entrainment; 
(□)  reactor  network  model  with  entrainment  estimated  from 
the  measured  data  (see  Fig.  2);  (+)  N02-reactor  network 
model  with  entrainment  estimated  from  the  measured  data 
(see  Fig.  2). 

above  analysis  and  run  5  temperature  plot,  shown  in 
Fig.  3  (left),  clearly  demonstrate  that  the  temperature 
increase  within  the  fuel  jet  results  almost  exclusively 
from  the  entrainment  of  the  hot  combustion  products. 
The  process  taking  place  within  the  fuel  jet  can  be 
perhaps  named  “preconditioning”  of  the  fuel.  This 
process  results  in  a  fuel  that  is  highly  diluted  with 
CO2  and  H2O  and  is  preheated  to  a  temperature  that 
is  typically  only  200-300  K  lower  than  the  temper¬ 
ature  of  around  1800  K  under  which  the  main  fuel 
oxidation  reactions  proceed  in  reactor  1 1.  In  run  5,  we 
have  calculated  that  such  a  preconditioned  fuel  is  of 
the  following  composition:  CH4  =  5.3%,  H2  =  0.5%, 
C02  =  13.5%,  H20  =  27.8%,  and  N2  =  50.9%.  The 
reactor  network  model  shows  no  NO  formation  within 
the  fuel  jet  and  its  presence  there  is  exclusively  due  to 
the  entrainment.  The  N02,  which  is  present  in  signif¬ 
icant  quantities  within  the  fuel  jet,  is  formed  from  the 


entrained  NO  via  reactions  N0  +  0  +  M  =  N02+M 
and  NO  +  H02  =  N02  +  OH  (the  first  reaction  pro¬ 
vides  around  two-thirds  of  the  N02  produced). 

6.  Conclusions 

We  have  performed  a  further  analysis  of  the  semi¬ 
industrial-scale  experiments  [1,7]  on  mild  (flameless) 
combustion  of  natural  gas.  The  experimental  burner 
has  featured  a  strong  (high-momentum)  oxidizer  jet 
and  two  weak  (low-momentum)  natural  gas  jets.  It  has 
been  observed  that  little  combustion  occurs  in  the  fuel 
jet.  The  temperature  rise  and  the  increasing  NO  con¬ 
centrations  along  the  fuel  jet  are  almost  exclusively 
due  to  the  entrainment  of  hot  combustion  products. 
It  has  been  proven  that  the  failure  of  a  number  of 
RANS  models  in  predicting  the  structure  of  the  fuel 
jet  is  a  result  of  error  in  predictions  of  the  entrain¬ 
ment  and  therefore  is  not  related  to  any  chemistry 
submodels.  The  novelty  of  this  paper  is  in  showing 
that  little  combustion  occurs  within  the  natural  gas 
jet.  The  temperature  increase  within  the  fuel  jet  re¬ 
sults  almost  exclusively  from  the  entrainment  of  the 
hot  combustion  products. 

The  process  taking  place  within  the  fuel  jet  can  be 
named  “preconditioning”  of  the  fuel  and  it  results  in 
a  mixture  containing  around  6%  of  combustibles,  2% 
of  oxygen,  and  92%  of  inert  gases  (C02,  H20,  N2). 
The  mixture  does  not  ignite  until  it  merges  into  the 
oxidizer  jet.  Consequently,  there  is  no  NO  formation 
within  the  fuel  jet  and  its  presence  there  is  exclusively 
due  to  the  entrainment. 
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